Indications of a shallow potential in ^^Ca+^^Zr fusion reactions 



H. Esbensen and C. L. Jiang^ 

^Physics Division, Argonne National Laboratory, Argonne, IlUnots 60439 
(Dated: June 12, 2009) 

Fusion data for ''**Ca+^^Zr are analyzed by coupled-channels calculations. Puzzling features of a 
previous analysis are eliminated by applying a potential that has a shallow pocket in the entrance 
channel. Thus the observed S factor for fusion, which develops a maximum at low energy, can be 
reproduced fairly well. The high-energy data can also be accounted for but that requires the use of 
a weak, short-ranged imaginary potential that absorbs the incoming flux near the location of the 
minimum of the potential pocket. Predictions of the fusion hindrance in other Ca-|-Zr systems are 
made and are compared with the systematics that has been developed previously. 
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INTRODUCTION 

In the field of heavy-ion fusion reactions two major is- 
sues have been discussed in recent years. One issue is 
how to explain the hindrance of fusion at low energies, 
that is the steep falloff of the measured fusion cross sec- 
tions, which has been observed in many heavy-ion sys- 
tems at energies far below the Coulomb barrier (see Ref. 
[H and references therein.) The hindrance is often so 
strong that the S factor for fusion develops a maximum 
at low energy, and it has been a challenge to reproduce 
this behavior in coupled-channels calculations. The other 
issue is the suppression of high-energy fusion data, which 
is observed when the data are compared to conventional 
one-dimensional or coupled-channels calculations ^] . 

A practical solution to the problems mentioned above 
has been to adjust the parameters of the Woods-Saxon 
shaped, ion-ion potential that is used in the calculations. 
Thus it has been realized that by using a large diffuseness 
(in combination with a deep potential and a small radius 
parameter) one can often resolve the problem of a sup- 
pression at high energy [1]. However, the diffuseness that 
is required to reproduce the high-energy data, for exam- 
ple, in the fusion of ^^0-|-^°*Pb , is different from the 
value that is needed to fit the low-energy data 0] . More- 
over, the values are often unrealistic and much larger 
than the empirical diffuseness which has been extracted 
from analyses of elastic scattering data or obtained from 
double- folding potentials Q • It has recently been shown 
that a large diffuseness is also inconsistent with quasi- 
elastic scattering data [1]. 

A good example on a system where both of the 
above mentioned problems exist is the fusion data for 
■^^Ca-l-^^Zr Q. The problems were solved in the origi- 
nal analysis of the data by using a large diffuseness, a 
= 0.85 fm, instead of the empirical value, a = 0.677 fm 
(Eq. (III.44) of Ref. fij). The measurements ^ were 
later supplemented with two new low-energy data points 
and they confirmed that the 5* factor for this system 
does develop a maximum at low energy. 

Instead of using an ion-ion potential with an anoma- 



lously large diffuseness, the double-folding technique is 
applied to construct an ion-ion potential which has a re- 
alistic behavior at large distances and is adjusted at small 
distances between the reacting nuclei to produce a shal- 
low pocket in the entrance channel potential Q. The 
potential used is the so-called M3Y-|-repulsion potential 
and the construction of it is described in Ref. Q. The 
basic feature of this potential is that the existence of a 
shallow pocket in the entrance channel makes it possible 
to explain the steep falloff of the fusion data at low en- 
ergies, and it also helps explaining the measured fusion 
cross sections at high energies [J S However, in 

order to reproduce the high-energy data, it is often nec- 
essary to include a short-ranged, imagi nary potential in 
the calculations of the cross sections [9j, |lO| . 

The M3Y-|-repulsion potential is applied in coupled- 
channels calculations of the fusion cross section for 
''^Ca-l-^^Zr. The purpose is to demonstrate that the en- 
ergy dependence of the measured S factor for fusion can 
be reproduced at low energies by employing this poten- 
tial. It is also shown that the calculated cross sections at 
high energies are sensitive to the nuclear potential and to 
the nuclear structure input, and it is possible to repro- 
duce the high-energy data by using a realistic structure 
input and employing the M3Y-|-repulsion potential, in 
combination with a short-ranged imaginary potential. 

Finally, the hindrance in the fusion of other Ca-|-Zr 
isotopes at low energies is discussed and compared with 
the systematics we have developed previously [Ij], and 
predictions are made of the energy where the S factor 
for fusion is expected to develop a maximum in those 
cases where it has not yet been observed experimentally. 

ION-ION POTENTIAL 

The construction of the double-folding ion-ion poten- 
tial for '^^Ca+^^Zr is based on the effective M3Y inter- 
action which is supplemented with a contact term for 
the exchange part [8|. The densities are parametrized as 
fermi functions with the radius 3.675 fm and diffuseness 
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FIG. 1; (color online) The M3Y+repulsion potential (red 
solid curve) is compared to the Woods-Saxon potential (mid- 
dle dashed curve) and to the pure M3Y-|-exchange potential 
(lowest black dashed curve). The diamonds show the poten- 
tial with large diffuseness that was used in Ref. 0]. The 
energy of the compound nucleus ^^''Nd and the energy band 
where the S factor for fusion has a maximum (« 88 MeV) are 
also indicated. 



TABLE I: Structure input for ''^Ca [Tl| and ^''Zr [H]. 



Nucleus 


A'' 


(MeV) 




(fm) 


(fm) 


*«Ca 


2+ 


3.832 


0.102 


0.126 


0.190 




2PH(2+) 


4.297 


0.082 


0.101 


0.101 




3" 


4.507 


0.203 


0.250 


0.190 


96 Zr 


2+ 


1.751 


0.079 


0.123 


0.123 




3" 


1.897 


0.295 


0.457 


0.457 



The entrance channel potential which was used in Ref. 
[sl (with diffuseness a — 0.85 fm) is indicated by the 
diamonds in Fig. [1] This potential is deeper than the 
other two potentials mentioned above but none of them 
are deeper than the ground state energy of the compound 
nucleus "'^'''^Nd, which is indicated by the lowest horizontal 
line. The broad horizontal band near 88 MeV indicates 
the energy where the experimental S factor for fusion 
develops a maximum. 

COUPLED-CHANNELS CALCULATIONS 



0.525 fm for ^Ca, and the radius 5.02 fm and diffuseness 
0.55 fm for ^^Zr. These two densities produce the RMS 
matter radii of 3.45 and 4.39 fm, respectively, which are 
consistent with the measured charge radii. 

The M3Y double-folding potential (including ex- 
change) is illustrated by the lowest dashed curve in Fig. 
[TJ It is unrealistic at short distances between the two 
nuclei because it is much deeper than the ground state 
energy of the compound nucleus. The M3Y interaction 
has therefore been corrected by supplementing it with 
a repulsive term as described in Ref. The densi- 

ties that are used in calculating the associated repulsive, 
double-fording potential have a sharper profile, with the 
diffuseness aicp=0.4 fm, but the radii are the same as 
used above to calculate the ordinary M3Y double-folding 
potential. The strength of the repulsive contact term, 
Vrep = 507 MeV, was adjusted to produce the nuclear 
incompressibility K = 223 MeV. 

The entrance channel potential for the M3Y-|-repulsion 
potential is illustrated by the upper solid curve in Fig. [T] 
It has a rather shallow pocket inside the Coulomb barrier. 
Also shown is the entrance channel potential obtained 
from a Woods-Saxon (WS) potential, which is similar to 
the empirical Akyiiz-Winther potential Q. The slightly 
adjusted parameters of this WS potential are Vo = -77.18 
MeV, R = 9.686 fm, and a = 0.677 fm. They produce the 
same Coulomb barrier height, Vcb = 95.4 MeV, as the 
M3Y-|-repulsion potential but the pocket is deeper. The 
fusion cross sections obtained in the coupled-channels cal- 
culations that are based on these two potentials are es- 
sentially the same at energies close to the Coulomb bar- 
rier but differ at energies far below the Coulomb barrier. 
This will be discussed in the next section. 



The nuclear structure input to the coupled-channels 
calculations associated with the excitation of the 2+ and 
3~ states in projectile and target is shown in Table I. 
The calculations include all onc-phonon, two-phonon, 
and mutual excitations of these onc-phonon excitations. 
That is a total of 15 channels in the rotating frame ap- 
proximation [8] . The coupling between the one- and two- 
phonon states is assumed to be vibrational (harmonic), 
except for the quadrupole excitation of ^*Ca for which we 
use a weaker one- to two-phonon coupling strength (see 
Table I) which is based on the known Oj 2^ transition 
strength [H]. 

The conventional assumption — f3'~^ does not 
work so well for '*^Ca. The empirical nuclear coupling 
strengths that will be used are shown in the last column 
of Table I. They are expressed in terms of the values of 

and were determined by analyzing measurements of the 
inelastic scattering of ^^O on ''^Ca 

The fusion cross sections are determined by imposing 
ingoing-wave boundary conditions at the position of the 
minimum of pocket in the entrance potential. The re- 
sults of coupled-channels calculations are compared in 
Fig. [2] to the data of Refs. [i,[3]- The top dashed (green) 
curve is based on the empirical Woods-Saxon potential 
discussed in the previous section. It is seen to reproduce 
the data quite well, except at the very lowest energy. The 
overall x^/iV (based on statistical errors only) is about 4 
for this calculation. 

The solid (red) curve in Fig. [2] is the coupled-channels 
result obtained with the M3Y-f repulsion potential. It 
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FIG. 2: (color online) Measured fusion cross sections 
for«Ca+s6Zr 0,0 are compared to coupled-channels calcu- 
lations based on the Woods-Saxon (WS, green dashed curve) 
and M3Y-|-repulsion (red solid curve) potentials. The low- 
est dashed (black) curve is the no-coupling limit for the WS 
potential. 
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FIG. 3: (color online) Barrier distributions obtained from the 
fusion cross sections shown in Fig. fusing AE = 2 MeV (open 
circles for AE — 1 MeV). The red solid curve was derived from 
the M3Y+repulsion calculation show in Fig. [2] The blue, 
dashed curve was obtained from a calculation that includes 
up to 5 phonons of the octupole excitation in ^^Zr. 



reproduces the data very well at low energies but is below 
the data at high energies. The latter problem can be 
solved as mentioned earlier by employing short-ranged 
imaginary potential. The explicit form used below is 



W{r) = 



-10 MeV 



1 -I- exp((r - rp)/aw) ' 



(2) 



where Vp is the radial distance between the reacting nuclei 
at the minimum of the pocket, and the diffuseness is 
set to 0.2 fm. The effect of the imaginary potential, Eq. 
dH), is illustrated below (see Fig. O) 



Barrier distribution 

The effect of couplings on the calculated fusion cross 
section is very strong. This can be seen by comparing 
the coupled-channels calculations to the no-coupling limit 
which is shown by the lower dashed curve in Fig. O The 
effect is to shift the calculated cross section by roughly 
4-5 MeV to lower energies. 

The couplings produce some interesting features in the 
barrier distribution for fusion which is defined as 
the second derivative of the energy- weighted fusion cross 
section with respect to the center-of-mass energy iJ, 



(3) 



The barrier distributions, which arc shown in Fig. [3l were 
calculated by the finite difference method using the same 
energy step AE = 2 MeV that was used in the earlier 
work 0. The measured distribution is broad and con- 
tains several peak structures. The red solid curve in Fig. 
[3] is the distribution derived from the coupled-channels 



calculations that are based on the M3Y-|-repulsion po- 
tential. It is broad and consists of two strong peaks but 
the shape is not in good agreement with the data. 

It was demonstrated in Ref. that some of these 
peak structures in the measured barrier distribution can 
be reproduced by including up to three-phonon excita- 
tions of the octupole mode in ^^Zr. The blue dashed 
curve in Fig. [3] is the barrier distribution one obtains by 
considering up to five phonon excitations of the octupole 
mode in ^^Zr in the coupled-channels calculations. The 
other excitations shown in Table I were not included in 
these calculations for practical reasons. The shape of the 
distribution (the blue dashed curve) seems to trace the 
four peak structures of the measured distribution rather 
well, and it appears to be in better agreement with the 
data than the three-phonon calculation presented in Ref. 
6] . It is obvious from this comparison that the calculated 
barrier distribution is very sensitive to multiphonon ex- 
citations of the soft octupole mode in ^^Zr. 

The fusion cross section obtained in the calculation 
that is based on up to five-phonon octupole excitations 
has an average deviation of about 10% from the data 
(i. e., we need a systematic error of 10% in order to ob- 
tain a psr point of 1), whereas the solid curve in Fig. 1 
deviates on average by 7%. It is possible that one would 
be able to improve the overall fit to the data by including 
more channels in the calculation. However, such a calcu- 
lation would probably not be very realistic because of the 
sensitivity to the empirically poorly known multiphonon 
excitations. The main point we want to make is that the 
standard calculation, which is based on couplings to one- 
and two-phonon excitations, provides a reasonable fit to 
the data and by applying the M3Y-|-repulsion potential 
it is possible to reproduce the fusion hindrance observed 
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FIG. 4: (color online) (a) The measured S factors for the 
fusion of ^®Ca+^''Zr ^3, i3] are compared to the two coupled- 
channels calculations shown in Fig. [2] The value of 770 was 
set to rio = 74.58. The blue dashed curve is the prediction 
of Eq. ([8| and the parameters of Table II. The logarithmic 
derivatives, Eq. are shown in (b). The top short-dashed 
curve is the constant S factor limit, Eq. ((6)1; the long-dashed 
curve is the no-coupling limit for the WS potential. 



at energies far below the Coulomb barrier. 



the fusion probability. 

The coupled-channels calculations that are based on 
the M3Y-|-repulsion potential provide the best descrip- 
tion of the low-energy S factor data shown in Fig. |4^. 
In fact, the S factor data have just reached a maximum 
at the lowest energies, which is a characteristic feature 
for the onset of the fusion hindrance phenomenon in sys- 
tems with a negative Q- value for fusion, and this feature 
is reproduced by using a shallow entrance channel poten- 
tial. The extrapolation to very low energies, which will 
be discussed later on, is also shown in the figure. It has 
a nearly parabolic shape and the lowest data points are 
seen to fall on the top of this curve. 

Another way to illustrate the fusion hindrance is to 
show the logarithmic derivative of the energy weighted 
fusion cross section. 



HE) 



dlniEaf) _ 1 d{Eaf) 



dE 



Ea 



dE 



(5) 



This quantity is shown in Fig. The calculations 

based on the Wood-Saxon potentials are seen to reach 
a local maximum at low energies. The data and the cal- 
culations based on the M3Y-|-repulsion potential, on the 
other hand, keep increasing steeply with decreasing en- 
ergy. They both intersect the upper 'constant S factor' 
line, which is the logarithmic derivative one obtains by 
assuming a constant S factor 1], 



0.495Z1Z2VM 



(6) 



Here ji = AiA2/{Ai + A2) is the reduced mass number. 
At the point where L{E) intersects with Lcs{E), the as- 
sociated S factor will have a maximum (see Fig. 2^.) 



Behavior at low energy 

The low-energy behavior of the fusion cross sections is 
illustrated in Fig. 2^ in terms of the S factor for fusion, 
which is defined by, 

SiE)^afiE) E exp(27rr^iE)y (4) 

where E is the center-of-mass energy, rj = ZiZ2e'^ / (hv) 
is the Sommerfeld parameter, and v is the relative veloc- 
ity. It is seen that the S factor for the coupled-channels 
calculations that are based on the Woods-Saxon poten- 
tial keeps increasing with decreasing energy. The cal- 
culation based on the M3Y-|-repulsion potential, on the 
other hand, develops a maximum before it goes to zero 
near 86.2 MeV, which is the energy of the pocket mini- 
mum for the M3Y-|-repulsion potential shown in Fig. [T] 
The calculated fusion cross section must go to zero at 
this energy because we use a real potential at low energy 
and use ingoing-wave boundary conditions to determine 



Behavior at high energy 

The behavior of the fusion cross sections at high en- 
ergies is illustrated in Fig. \5[ The lowest dashed curve 
is the result of the coupled-channels calculation that is 
based on the M3Y-|-repulsion potential. This calculation 
is suppressed compared to the data, as already pointed 
out in the discussion of Fig. [21 However, the discrepancy 
with the data can essentially be eliminated by includ- 
ing in the calculation the weak, short-ranged imaginary 
potential, Eq. [51 This is illustrated by the red solid 
curve. The result obtained with the Woods-Saxon po- 
tential (with diffuseness a — 0.677 fm) and the structure 
input from Table I is shown by the blue dashed curve 
which is difficult to see because it is essentially identical 
to the solid curve in Fig. [51 This remarkable feature is 
the reason why the two calculations give essentially the 
same quality fit to the data as mentioned earlier. 

Good agreement with the high energy data could only 
be achieved in Ref. Q by using a large diffuseness in the 
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FIG. 5: (color online) Linear plot of the fusion cross sections 
for **Ca+®®Zr @|. The upper dashed curve is based on the 
WS potential with diffuseness a — 0.677 fm and the structure 
input from Ref. The next (blue) dashed curve is based 
on the same potential but stucture input is from Table I. The 
red solid and the lower dot-dashed curves are the coupled- 
channels calculations, based on the M3Y-|-repulsion potential 
and the structure input of Table I, with and without the effect 
of the short-ranged imaginary potential iW{r), respectively. 



Woods-Saxon potential. This seems to contradict the re- 
sults discussed above, where the Woods-Saxon potential 
with 'normal' diffuseness provided an excellent descrip- 
tion of the high-energy data when applied in coupled- 
channels calculations. The reason for this apparent con- 
tradiction is that the structure input of Table I is different 
from the structure input used in Ref. in particular 
for ''^Ca. The coupled-channels calculations were there- 
fore repeated using the structure input from Ref. Q and 
a Woods-Saxon potential with the 'normal' diffuseness a 
= 0.677 fm. The result is shown by the upper dashed 
curve in Fig. [S] It is seen that the data are suppressed 
compared to this calculation, in accord with the findings 
in Ref. 



FUSION HINDRANCE IN CA+ZR SYSTEMS 

Many fusion reactions involving medium heavy nuclei 
have been measured at Legnaro. Here we focus on the 
measurements that have been performed with different 
isotopes of calcium and zirconium The pur- 

pose is to extract information about the fusion hindrance 
phenomenon at extreme subbarrier energies and to com- 
pare with the results that have been obtained for other 
systems 

The fusion hindrance for systems with negative Q- 
value can be characterized by the energy Eg where the 
S factor for fusion has a maximum at low energy 
Although many measurements do not reach such low en- 
ergies, it is still possible to determine Es by extrapo- 
lation of the existing data. It was found that the log- 



TABLE II: Parameters for the fusion hindrance behavior, 
Es, La, Ao, Bo and as, obtained for the Ca-|-Zr systems. The 
last column is the reference energy 

defined in Eq. 



System 


Es 


Ls 


Ao 


Bo 


as 


EV' 




MeV MeV"^ 


MeV"^ 


MeV^/^ 


mb 


MeV 


«Ca+^°Zr 


93.2 


2.31 


-17.1 


4222 


0.37 


92.8 


«Ca-H9^Zr 


86.4 


2.61 


-11.6 


3585 


3.8E-3 


93.2 


«Ca+''«Zr 


83.1 


2.77 


-7.43 


2783 


3.7E-5 


93.4 


^Ca-H^OZr 


90.0 


2.59 


-12.9 


3775 


0.025 


96.7 


^Ca+^'fiZr 


88.1 


2.71 


-13.8 


4528 


0.030 


97.4 



arithmic derivative of the energy-weighted fusion cross 
section, the L{E) defined in Eq. ([5|), has a nearly lin- 
ear dependence on energy at energies around the S fac- 
tor maximum P, [l^. A slightly different (non- linear) 
parametrization is used in the following [l^ , namely. 



L{E) = Ao 



Bn 



[E + Q)^p 



(7) 



where Q is the Q value for the fusion reaction, and 
we assume that Np=3/2. This expressions implies that 
L{E) — > oo for E — > —Q in agreement with the constraint 
that the cross section must go to zero when the center-of- 
mass energy E approaches t he g round state energy, —Q, 
of the compound nucleus [H, |l8| 

Once the parameters Aq and Bq in Eq. [7] have been 
determined from the low energy data, it is straight- 
forward to extrapolate and find the energy Es where 
L{E) intersects the logarithmic derivative Lcs{E) for a 
constant S factor, Eq. ([6]). Results for the systems 
40Ca-|-9°'94'9^Zr and ^^Ca-t-^^Zr using this extrapolation 
method are shown in Fig. [6^. The parameters obtained 
by fitting the low-energy L{E) data and the extrapolated 
values of Eg are collected in Table II. There we also give 
the value which is the logarithmic derivative at the 
energy E^. 

The analytic expression for the cross section one ob- 
tains from the parametrization ([7]) of L{E) is. 



a{E) 



Bo 



El 

E 



e^]i{An{E-Es 



{Np-l){Es + Q) 



, Es + Afp-i _ 
^E + Q> 



(8) 



The only unknown parameter in this expression is the 
normalization factor Us which is the fusion cross section 
at the energy E^. It can be obtained by normalizing to 
the low-energy fusion data and the resulting values are 
given in Table II. The S factors that have been deter- 
mined in this way are shown by the curves in Fig. [6Jd. 
They have a characteristic, nearly parabolic behavior at 
low energy. The measured S factors have not reached 
a maximum for most of these systems but the data for 



6 



3 ^ 



> 9 





102 

CN 
"ft 

g 10-2 

w 
e 



10- 



I I I I I I I I I I I I I I I I I I 

(a) 40ca + 90zr 




I I I I I I I I I I I I I 



40Ca + 96zr\ 




80 84 88 92 

E - Eg (MeV) 



FIG. 6: (color online) Logarithmic derivative (a) and S factors 
(b) for the fusion of «Ca+''°'^'*'^SZr [l5, 16] and ''«Ca+*Zr 
The curves in (a) are fits to the measured L{E) using 
the parameterization Eq. ((Tj. The upper, nearly horizontal 
dashed lines show the constant S factor limit, Eq. ([6]). The 
curves in (b) are the S factors predicted by Eq. (|8]l and the 
parameters in Table II. The parameters Eq and ryo used in 
the plot are Eo = 0, 1, 3, MeV and 770 = 68.59, 70.87, 
71.71, 73.98 for the systems '*''Ca+'"''^'''^^Zr and ''*Ca+^°Zr, 
respectively. 



^"Ca+^^Zr are close (although the lowest point has a 
somewhat unusual behavior.) The S factor for the sys- 
tem ^^Ca+^^Zr, which was already shown in Fig. 2^, did 
reach a maximum at the lowest energy, and the extrapo- 
lation obtained from the parameters in Table II was also 
shown there. 

We have previously studied the systematics of the en- 
ergy Eg for other heavy-ion systems (Refs. [J, [l3|) and 
found that the value of L^, i. e., the logarithmic deriva- 
tive of the energy weighted cross section at the energy 
Eg, has a nearly constant value of L™^ = 2.33 MeV~^ 
for systems that consist of closed shell nuclei or of nuclei 
that are relatively stiff. Inserting such a constant value 
of Lg and the energy Eg into Eq. ([5]), it is seen that 
the system dependence of Eg is governed by the entrance 
channel parameter 



C = ZiZ2y/JI = Z,Z2VAlA2/{Ai+A2 




^Ca+SfiZr" 



4100 4200 4300 ,„ 4400 4500 4600 

FIG. 7; The energy Es for Ca+Zr and Ni+Ni systems. The 
solid curve is the reference value, Eq. ([9]). The open symbols 
are extrapolated values. The solid symbols are for systems 
where a maximum in the S factor has been observed. 



The expression one obtains for Eg is the reference energy 

Ef = {0A95C/Lf)^/^ (MeV), (9) 

which is illustrated by the solid curve in Fig. [7l For 
systems consisting of soft or open shell nuclei, the value 
of Eg is usually lower than the reference value E^^^ , The 
deviation of Eg from El°^ is a measure of the softness the 
systems [Uliili^. 

The values of Eg extracted from the data for the dif- 
ferent Ca-I-Zr systems [1, 0, El, are shown in Fig. 
[71 Also shown are the values obtained previously for the 
Ni-|-Ni systems, consisting of different nickel isotopes. It 
is seen that the extrapolated value of Eg for ^°Ca-|-^°Zr 
falls on the solid line, consistent with the closed shell na- 
ture of these two nuclei. The results for the two soft and 
open shell systems, "'^Ca-l-^^Zr and ^'*Ni-|-^'*Ni, are simi- 
lar and fall far below the solid curve. The value of Eg for 
^*Ca-|-^^Zr, for example, is 9.3 MeV below the reference 
value i?"^, according to Table II. 

The values of Eg are also far below the reference value 
£:f f for some of the other Ca-|-Zr systems. This is partly 
caused by the softness of the reacting nuclei or by transfer 
reactions. Stefanini et al. in Ref. [16| pointed out that 
neutron transfer with positive Q value should affect the 
fusion of ^°Ca-f ^''Zr and ''°Ca-F^^Zr and enhance it below 
the Coulomb barrier. Consequently, the fusion hindrance 
would be pushed to lower energies. 

There is a correlation between the ratio Eg/El°^ and 
the number of valence neutrons in the two reacting nuclei, 
which was first pointed out in Ref. [20j for nickel induced 
fusion reactions. A similar trend is observed in Fig. [7| 
for the Ca-|-Zr systems. In order to see the interplay 
between neutron transfer with positive Q value and the 
fusion hindrance behavior, it would be very interesting to 
measure the fusion of ^"Ca-l-^'^Zr and ^"Ca-l-^^Zr down to 
much smaller cross sections. 
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CONCLUSIONS 

We conclude that there are strong indications of a fu- 
sion hindrance at low energy in the data for '*^Ca+^®Zr. 
The S factor does reach a maximum at the lowest ener- 
gies and there is a preference among the coupled-channels 
calculations presented here for those that are based on a 
shallow potential in the entrance channel. In this re- 
spect, the fusion data exhibit the same characteristics 
that have been observed in many other heavy-ion sys- 
tems 

The systematics of the fusion hindrance in other 
Ca-|-Zr systems, where it has not yet been confirmed ex- 
perimentally, was discussed and predictions were made of 
the energy Eg where the S factor for fusion is expected to 
develop a maximum. This energy can be far below the 
energy predicted by the empirical formula, which was 
developed for systems consisting of stiff or closed shell 
nuclei. The deviation of Eg from the empirical formula 
seems to correlate with the number of valence neutrons 
in the reacting nuclei, which may reflect an increased 
softness or a stronger effect of neutron transfer. 

The behavior of the high-energy fusion data can also 
be reproduced by applying the shallow potential in the 
coupled-channels calculations but that requires an addi- 
tional weak absorption that acts near the minimum of 
the potential pocket. The calculated, high-energy cross 
sections are sensitive to the nuclear structure input pa- 
rameters. By applying the structure input which has 
been extracted from inelastic scattering data it is possible 
to reproduce the high-energy data, even with a Woods- 
Saxon potential that has a normal diffuseness. 

The energy dependence of the fusion cross section in 
the barrier region was also investigated. This region 
is best studied in terms of the extracted barrier distri- 
bution, which is usually very sensitive to the low-lying 
nuclear structure input in the coupled-channels calcula- 
tions. The calculated barrier distribution for ''^Ca-f^^Zr 
turned out to be very sensitive to multi-phonon excita- 
tions, in particular of the soft octupole node in ^^Zr. The 
best agreement with the measured barrier distribution 
was achieved by including excitations up to five phonons 
of this excitation mode. Unfortunately, the empirical 
knowledge of the nuclear structure at such high excita- 



tion energies is poor and many other reaction channels 
may also affect the barrier distribution. It is therefore 
a serious challenge to theory to determine a credible nu- 
clear structure input and explain the measured barrier 
distribution in detail. 
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